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Abstract 

The Kontsevich-Penner model, an Airy matrix model with a logarithmic po- 
tential, may be derived from a simple Gaussian two-matrix model through a 
duality. In this dual version the Fourier transforms of the n-point correlation 
functions can be computed in closed form. Using Virasoro constraints, we 
find that in addition to the parameters t n , which appears in the KdV hierar- 
chies, one needs to introduce here half-integer indices t n /2 ■ The free energy 
as a function of those parameters may be obtained from these Virasoro con- 
straints. The large N limit follows from the solution to an integral equation. 
This leads to explicit computations for a number of topological invariants. 
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1 Introduction 



In some recent articles [TJ [2j [3j SJ [5] , we have discussed a relationship between 
the Airy matrix model (Kontsevich model) and a Gaussian random matrix 
theory with an external source. The free energy for higher Airy matrix mod- 
els (degree more than three) is the generating function for the intersection 
numbers of the moduli space of p-spin curves [HI [7J IE]- We have shown that 
the Fourier transform of the n-point correlation function of the Gaussian 
random matrix is also a generating function of these intersection numbers 
with n- marked points [H H] . The reason for this remarkable agreement is a 
duality for expectation values of characteristic polynomials [2J. 

It is interesting to extend this duality to the case of the c = 1 [H] matrix 
model, i.e. models dealing with a ID chain of coupled matrices, and to 
understand the meaning of the topological invariants. In [3], we have shown 
that the time dependent random matrix theory, is equivalent to a two-matrix 
model. These matrices are N x N and one considers correlation functions 
involving k\ points for the first matrix, and k 2 for the second. After duality, 
we have shown that the edge scaling limit is an Airy matrix model with a 
logarithmic potential 



where B is a ki x k\ Hermitian matrix, and k% is the number of characteristic 
polynomials in the two-matrix model [3]. In section 2, we recall the derivation 



This Airy matrix model with a logarithmic potential, the so called Kontsevich- 
Penner model, and its generalization to higher Airy matrix model (degree 
larger than three) have been discussed in the literature. For the general 
(p + l)-valent vertex model V(B) = B p+1 / (p + 1), the partition function 



was considered by Mironov et al [TO], in particular for small A, through 
character expansions. The parameter A separates (i) the strong coupling 
and (ii) the weak coupling regions. The transition between (i) and (ii) is 
similar to the Brezin-Gross-Witten transition [111 112] in the unitary matrix 
model. Indeed, it can be shown that the partition function Z in (II. 2ft with 
p = —2 is equivalent to the unitary matrix model (see Appendix B). For 
p = 1, it has been studied in [131 El- 

The model (ll.2p has been considered with polynomial vertices 




(1.1) 



of (fLTj). 




(1.2) 



V(B) = y £t n B n . 



(1.3) 
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as a generating function of tachyon amplitudes in the c=l string theory [TSJ 
[Tof [T7] . The existence of a two dimensional black hole [TTJJ [201 HD] has been 
discussed in related matrix models. We also note that some time evolution 
problems, such as crystal growth or non-intersecting random walks, share 
interesting universal features described by the edge singularity of a random 
matrix theory [211 [221 1231 122 [221 12E] ■ 

In this paper, in view of these interesting applications, we study the Airy 
matrix model, and higher Airy matrix models, with a logarithmic potential, 
in detail. We first obtain the Virasoro constraints for the partition func- 
tion Z of (11. ip ; this gives a series expansion in terms of the parameters 
t n = trA~ n ~2 ; constructed from the source matrix A. The parameters t n 
characterize the intersection numbers of the moduli space of curves and the 
KdV hierarchies, following Witten's well-known conjecture [7J [8]. The Vi- 
rasoro constraints are derived from the equations of motion for this matrix 
model. The remarkable difference with the Kontsevich model [30J is that 
the equation of motion becomes here a third order differential equation, due 
to the presence of the logarithmic term. This differential equation leads to 
the appearance of a new series of parameters t n / 2 (n/2 is a half- integer) in 
addition to the t n (n is an integer) of the KdV hierarchies. It is interest- 
ing to analyze the role of these new variables t n /2, absent from the usual 
KdV hierarchies. These parameters t n /2 (n/2 half- integer) correspond to the 
Ramond sector in string theory, decoupled from the Neveu-Schwarz sectors 
in the Kontsevich model or in the theory of intersection numbers of Rie- 
mann surfaces. Parameters similar to t n /2 appear also in the antisymmetric 
Gaussian random matrix theory [27J- 




In the large N limit the equation of motion leads to a Riemann-Hilbert 
integral equation. We have verified the consistency of the solution with the 
results derived from the Virasoro constraints. The free energy is expressed 
through the parameters t n and is. Starting from the duality relation between 
the Gaussian two-matrix and the Kontsevich-Penner model, we consider the 
Fourier transform of the correlation functions in the two-matrix model. The 
terms i« are generated by the correlations between the edge of the spectrum 
and the bulk. The edge behavior and the bulk behavior obey different scal- 
ing laws. This paper is organized as follows; in section 2, we recall how a 
logarithmic potential is generated from a duality relation for characteristic 
polynomials in the Gaussian two-matrix model. In section 3, the equations 
of motion are derived for the Kontsevich model (p=2) with a logarithmic 
potential. The Virasoro constraints are obtained, as differential equations 
for to, ti and t\. From these differential equations, we construct the series 
expansion of the free energy F in terms of the t n and is. In section 4, 
the integral equation (Riemann-Hilbert problem) is solved for the Airy ma- 
trix model with a logarithmic potential. In section 5, we discuss the replica 
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method which provides explicit results for one marked point. The corre- 
sponding correlation function of the two-matrix model is discussed and the 
result of its Fourier transform is compared with the result of the free energy 
obtained from Virasoro constraints. The section 6 is devoted to discussions. 
The formula which expresses the p-th derivatives with respect to the external 
source matrix A in terms of its eigenvalues is presented in the appendix A. 
In Appendix B, the equivalence between the unitary matrix model with an 
external source and (p = — 2)-higher Airy matrix model with a logarithmic 
potential is briefly sketched. 

2 A logarithmic potential 

For the one-matrix model, we have shown that the Kontsevich model is 
related to a matrix model at the edge of its spectrum, through a duality 
relation and the replica method. For the mathematical definition of the 
intersection numbers of the moduli space of curves, we refer to [8] . It involves 
an integration over the compactfied moduli space M g>n with genus g and n- 
marked points, 

<r dl ---r dn >= I c 1 (£ 1 ) dl ---c 1 (C n ) d ". (2.1) 

where c\ is the first Chern class and £i is a cotangent line bundle at the 
i-th marked point. This definition of the intersection numbers has been 
generalized [6] to the moduli space of p-spin curves. The intersection numbers 
have now an additional spin- index like in r n j, in which j takes values from 
to p — 1. They are defined by 

<r dl , n ---r dn , n >= - f_ c D (V)c 1 (£ 1 ) d i---c 1 (£ n ) d ". (2.2) 

where cd(V) is a D-dimensional top Chern class of the vector bundle V, 
V = H 1 ^, T). There is a cover of Riemann surface, and the line bundle £ 
has p roots, 

£ ~ T 0P (2.3) 

where T is isomorphism class. There are (p — 1) roots corresponding to the 
Neveu-Schwarz sector and one to the Ramond sector in string theory. 

These intersection numbers with spin indices may be computed from the 
p-th higher Airy matrix model. Non-vanishing intersection numbers satisfy 
the condition, 

n 

Y,di + D = 3g -3 + n (2.4) 

i=l 
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In the limit p — > — 1 the top Chern class becomes the Euler characteristics 
with marked points. In this case, the higher Airy matrix model reduces 
to the Penner model (logarithmic potential), from which one computes the 
Euler characteristics Xg,n [311 [38] . In this section, we will discuss the duality 
between the Kontsevich model with a logarithmic potential (II -ip and the 
Gaussian two-matrix model with an external source. 

In previous papers jU [2H], we have discussed the time dependent (c = 1) 
Gaussian random matrix theory. We have shown that the correlation function 
at two different times is equivalent to the correlation function of a two-matrix 
model with the following probability distribution, 

P(M U M 2 ) = — exp( - -trMi - -trMf - ctrM 1 M 2 + trMiAi + tiM 2 A 2 \ 
Z V 2 2 / 

(2.5) 

where Z is a normalization constant. 

We consider the correlation function i 7 fe lj fe 2 of characteristic polynomials, 
which is defined by the (k\ + fc2)-point correlation function, 

Fk lM =< II det ( A « - M i) II det (/^ - M 2) > ( 2 - 6 ) 

a=l 0=1 

where the average < ■ ■ • > is computed with the distribution P in (12. 5p . 

Let us briefly review how this correlation function reduces to the Kont- 
sevich model with a logarithmic potential pE]. This Fki,k 2 is expressed as a 
Grassmann integral over i[) a (a = 1, ki), Xp(P = 1) ^2), 

F klM =< f d^dHxdx^ a{Xa ~ Ml)i>a+M ^' M2)xp > (2.7) 



The integration over the matrices Mi and M 2 in then easy with the Gaussian 
distribution (12. 5p . It yields quartic terms in ip and x- These terms may be 
expressed as integrals over Hermitian matrices Bi, B 2 and complex matrices 
D, . {ki x ki for the auxiliary matrix B\, k 2 x k 2 for B 2 , and k\ x k 2 for 
the complex matrix D). 

N - - r ( N iN - \ 

ex P[~ 2 (i - C 2) ^^J = J dB i ex Py - y tr5 i + ^ 1 _ c2 B ^j 

ex P[~ 2 (i N _^ XXXX\ = J dB 2 exp(^ - ytrS| + -^==B 2 xxj 

Nc — f / N fc - 

exp[ ^XXip] = / dDdD^expi - NtrD^D + v tr(D4>x + P>^X^) 

1 — c 2 J \ v 1 — c 2 

(21 

Note that we have traded the integrations over N x N matrices by integrals 
over matrices whose sizes are given by k\ and k 2 . One can then integrate out 
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the Grassmann variables ip and x, and F kljk2 is expressed as 

F klM = J rffi 1 rf J B 2 rf J Drf J D t e-T to (^+B 2 2 +2DtD )+ ivt r iog(i.x) (2 g) 



where 



with 



X = "i- c "2 ai-caa _ (2.10) 



a2— cai (i2~ ca\ 



(-Bl)aa' = (-Bl)aa' ~ Wl ~ C 2 \ a 5 aa > 



(2.11) 



We have assumed that the external source matrices A 1 and A 2 are multiple 
of the identity A ± — a± • I, A 2 = a 2 • I- Introducing the diagonal matrices Ai 
and A 2 

Ai = diag(Ai,...,A fcl ), A 2 = diag(//i, (2.12) 

we obtain 

F fcl , fc2 = e f (l- C 2 )tr(A 2 +Ai) J dB x dB 2 dDdrte~^ tr(B 1 2 + B|+2DtD)+iVt r lo g (l-X) 
x e -jAf v / l-c 2 tr_BiAi-jAfv / l-c 2 tr_B 2 A2 ^2 13) 

We now restrict ourselves to a 2 = 0. From the expression of X, we have 

9 (1-c 2 ) ~ 9 1-c 2 ~, 2(1 -c 2 ) t 
trX 2 = ^trR 2 ^tr£ 2 - 1 _ k rD^D 2.14 

If we take a 2 = 1 — c 2 , the quadratic term in B\ is cancelled. The quadratic 
term in B 2 does not vanish, it becomes — y(l — c~ 2 )tiB 2 . The quadratic 
term trD^D is also cancelled. 

We now denote 5 2 by B 1: B 2 . We find 

v3 i(l -c 2 ) 3/2 n3 3z(l -c 2 ) 3 / 2 „ , 3i(l -c 2 ) 3 / 2 . „ 
trX 3 = — i =^ — trB 3 i — tiDD^Bx + — ^ ^ — trD^DB 2 



(i\ ' " ca\ 



( /J «bi am 

c Hi 



Given the factor N in the exponent, the edge scaling limit under consid- 
eration corresponds to 

B 1 ~ 0(iV _ 5), B 2 ~ 0(N~*), D ~ O(aH) (2.16) 
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in the large N limit, since the quadratic term in B 2 does not vanish. In this 
limit most terms disappear ; for instance 

Ntv(D^DB 2 ) ~ N-t, (2.17) 

is negligible. Then, in the large N limit f 1 2 . 1 6 j) . after dropping the negligible 
terms, we obtain 

F klM = J dBl dB 2 dD^dDe~ imTB ^- imTB ^^^ 

(2.18) 

Since the matrix B 2 is decoupled, we can integrate it out. Then, dropping 
the contribution from the integral over B 2 , we find 

F klM = J dB 1 dD^dDe- itlB ^ tlB " +itWDfB ' (2.19) 

where we have absorbed the powers of N in a rescaling. 

We may now integrate out the matrices D and (D is a k\ x k 2 complex 
matrix); this yields a one matrix integral with a logarithmic potential, 

F klM = J dB ie ^ lB "- k2tTlogBl - itlBlAl . (2.20) 

where B\ is a k\ x k\ Hermitian matrix. If we chose the replacement B\ — > 
—iB, we obtain the model ( II. ip . 

3 Virasoro constraints 

The Kontsevich model with a logarithmic potential, 

Z = J dBe tr(-^ + AB + klo g B)^ (31) 

in which B is an Hermitian P x P matrix (we have replaced k\ by P and k 2 
by k). satisfies the trivial equations of motion, 

f dB J_ e ^-\B* + AB + k\o g B) = g (32) 

J oB ba 
from which one obtains readily 

(-dv)l + ( AT ^) oi + ^ + ^) z -°- (M) 
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Since Z is a function of the eigenvalues \ of A, one can trade this for 
differential equations in terms of these eigenvalues (see appendix A), 

d 3 Z v 1 _g d d d 

dX c 3 + ^ c A c - A/ dX c dX d K dX c + dX d } 

1 d d 1 d d 

'^(K-Xd^dXc ~ d\d Z + 2 £ c (A c -A e )(A e -A,) ( 9A; " dX~J 
8Z 

-A c — -(P + k)Z = (3.4) 

The zero-th order contribution for large A, is obtained from the shift 
B — > B + A 2 ; then keeping only the terms which grow for large A one finds 



/ 



dBe^ 11 +if trA5 +f trl °s A 
1 3 



e-^UXf (3.5) 



In the limit A, — > oo, the partition function reduces to Zq. Therefore, the 
partition function Z may be expressed as 

Z = Z g(X) (3.6) 

where g has an expansion in inverse powers of y/X : 

g = l + 0(X) (3.7) 
A2 

The Virasoro constraints (13. 4 ft lead to a sequence of equations, which fix 

the coefficients of the terms A c 2 . Let us thus write the Virasoro constraints 
in terms of the function g of (13. 61) . For this purpose, we have to substitute Zq 
into (13. 4p . The resulting equations are cumbersome. To avoid complicated 
and long expressions, we take the simple case of P = 2. Although this simple 
case is manifestly not sufficient to determine the expansion in terms of the 
t n , it is instructive and useful also for arbitrary P as shown below. Then for 
P = 2, the equations (13 .4p become 

<9 3 I , d d s/n d d x 
+ \ H7JT--7rH( 2 7JT- + 



d\\ Ai — A 2 d\i dX 2 dXi 8X 2 

1 d d d 

A 17 ^-(2 + A;))Z = (3. 



(X 1 -X 2 ) 2 "dX 1 dX 2 J dXi 

Using Z = Z g, we obtain the equations for g, 

dg dg d 2 g d 2 g d 2 g d 3 

a x g + a 2 — + a 3 — + + a 5 ^- + a 6 ^ + a 7 ^ 
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where 



ai 



l-2fc\ 1 /l — 2fc (2/c — l)(2/c — 5) 



+ 



V A1A2 

1 /5 - 24/fc + 12/fc 2 k 
+ "( 16 + 



Ai V 2yA2 
1\ 1 



2AJ 



16A^ 



(Jfe-l)(2Jfe-3) 



+ 



1 (k - 1 



+ 



Ai 2 V A 2 2 
-1 + 16A; - 12fc 2 



2" ' ^1 

J_f(2k- l)(2fc — 5)(2/c — 9) 

a? 



4A, 



64 



1 



32Ai 2 



/AT+ VA^ 



(3.10) 



1 

«2 = - 



1 - 2A; 
+ ^ ; + 



2^+^ 2 a/Ai - V^2 2A X A 2 A^a/A 



3A; 



+ 



+ 



4A12VAT+VA2 A X 2 VAT-VA2 4A 1 (v / AT-VA2) 



-(3.11) 



a 3 = 2Ai + 



3A; 



+ 



15 - 36A; + 12A: 2 1 - 2k 
~ ~ + 4A^A7 



16A^ 
3A; 1 



4A 1 2v / AT+v / A 2 \ X 2 VAT-VA2 4A 1 ( V / AT-VA2) 



04 = 3y Ai — 



3(1 -2k) 
4A^ 



+ 



j(3-12) 

(3.13) 

(3.14) 

(3.15) 

We now return to general P (not simply P = 2) and define the parameters 

1 



2VWAi + VA^) vAi(vAi - VA2) 
1 

a 5 = a 6 = 



Ai — A2 



a 7 



t. 



E 



1 Aj 



(3.16) 



in which n takes both integer and half-integer values ( n = 0, |, 1, |, • • • ). 
Note that only integers appear in the Kontsevich model. The appearance 
of half-integers is a characteristic of the present model with a logarithmic 
potential. The derivatives with respect to Xj are replaced by 



d 



0t„ d 1 ] 



9 



<9Aj <9t n 



2 y A/ + 5 dt r . 



(3.17) 
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y QMj)(n + |) _g vv (n + |)(m + |) 

aA, 2 ^ v+§ a*» r£ a,™ 3 <% n 0* m 

d 3 y^(n + 5)(n+§)(n+§) 9 



0A x d V A^* <%n 

(n+ |)(m+ ±)(2n + m + f) 2 



EE 



(n + |)( m +|)(j- + i) 03 

r ™ ^ A r+ m+ ^ 0^0^ 1 ' ; 



where n, m, j = 0, \, 1, |, 2, • • •. 

Returning now to P = 2, at lowest order in the 1/a/Ai expansion, a\ 
becomes 

1/1/1 1 „ LI 1 , \ , 1 , 

°'~7xr(4^ + 735 ) -2^ + ^) +0( aI ) (3 - 21) 



This may be expressed in terms of the t n in (I3.16P as, 



«i ~ 7/^4*0 - (3-22) 



9 1 ,1 , . , 

tt2 ^~-7S7 ( 2^ + i la7 + "' ) (3 ' 23) 

2 

<9 9 19 , 

a 3 ^- ~ 2A X — =— 3.24 

0Ai 0Ai vAi 0t o 

We have considered only the case P = 2, but this simple calculation is enough 

to determine correctly the coefficients of £q and ti, which are defined as the 

sum of A" up to i = P as f)3.16p . 

The coefficients do not appear yet at this order since the 

multiplications of derivatives in f !3.9p give higher orders in A^f 1 . Then, we 

obtain the first equation of order X t , 

- ? (n+ 5 )t » + 'l:) 9 = (3 - 25) 

n=0, 5 ,l,... 

Using F = logg, it becomes 

dF 1 „ fc „ . 1 0F . 

n=0,i,lr- 



For the next order Xj 1 , we need to evaluate a\ for N = 3, since the N = 2 
results are not sufficient to determine the coefficients of t n which appear in 
the equation. We obtain 

1 (A k M l 1 



y/Xx V 4 2 A2 A3 2v^A2A3, 
1 (\ ,., 1 1 , (-5 + 12* - 4fc 2 ) , 1 1 , 

/ 1 fcv/ 1 1 1 A , 

+ (-o + o)(TV + ^ + ~ )) ( 3 - 2? ) 

2 2 A 2 A 3 AfVA^ Afv^ 7 

From fl3.lip . f13.12p and (13.13p . we obtain the second equation, involving now 
a derivative with respect to the ti, 

„ d , k 9 k 1 9 f 1 , „ 

2- ton + -ti + ~* *1 t Hi ta kHz 

dti 4 2 2 u 4 5 16 2 4 2 

2 

- E ( 2 " + 1 )W|«r + ^ E ( n + 2 )W2 9T 

n=0±l,.. ra n=0±l,.. n 



2 ■ 

,2 



(3.28) 



2 . — :; Z t/Cjfc Z — 8 Z Z CTCiCTCjfc 

— I— J+fc=— f -t+7+fc=-§ 



The next order is proportional to Ai 2 , and we obtain 

„ d 1 3, 2 k- 1 
3 S7 - 16 - 4* + — *-** 

These equations determine the free energy F = logg as up to order 

o(x-l), 



1 3 1 1, 1,9 

F = — 1„ + — *i + -ktoti + -k 2 h 
12 48 2 5 4 



+ — + (— + —k 2 )t 2 + -fain* 1*1 + — to 3 i + (— + -k 2 )t t 2 + -ton*s 
24 1 V 192 16 ;i 4 2 1 24 I V 32 8 ; 4 °5 

+ -kHits + -(A; + A; 3 )t 5 
4 2 2 6 2 

+ — *n*2 + -ton*s + — tit 2 ! + (— + — k 2 )tlh 
64 6 ° 2 48 1 v 128 32 !Qi 

3 11 1 

+ T77to^it 2 + -MltiU + -(A; + fc 3 )t *Z + ^k 2 t tit 5 _ 
16 2 4 2 2 2 12 22 
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+ + a k2 )W2 + -k\t\ + -kt t\n + -kt Q titl 
32 8 4 282282 

1 73 43 3 1 1 

+ — ( — k H k 3 )tit 3 H k 2 t\t 2 + -kHUxtz + -(fc + fc 3 )t 1 * 5 

14 v 32 8 ; 2 3 ^i6 2 2 ^4 2 X 2^6 V ; 2 

,1 1,.., A, 1,,, 1, 105 607, 2 169, 4N 

+ + — fc 2 t? + (-fc + -k 3 )t 3 t 2 + - + k 2 + k A )t 4 

V 576 48 ; 1 V 8 4 ; 2 2 ^ 9 V 1024 128 64 ; 

+ -^tiiti (3.30) 
24 2 

This expression is consistent with the previous result [4j. Note that the 

parameters t n with the half-integers n, ti.ta,... appear together with the 

2 2 

coefficients proportional to k. When k goes to zero, the free energy F of 
(I3.30P reduces to the Kontsevich free energy. Another remarkable propertiy 
of f)3.30p is that when k is of order P, many terms are of the same order in the 
large P limit. The leading order is P 2 which gives genus zero contributions. 
We will discuss the large P limit in a later section from a different approach 
based on integral equations. 

To express these equations in compact form, it is convenient to introduce 
the differential operators obtained as follows 



J m {1) (x) = jL- mx _ m , (m = ..,-2,-1,0,1,2,...) (3.31) 

ox m 

and x m = for x > 0. We define (k > 1) from J$ as 

4?= E ■■J < i 1) jf ) -- (3-32) 

i+j=m 

where : • • • : means normal ordering, i.e. pulling the differential operator to 
the right. Then we obtain 

d 2 
Jj 2) = E + 2 E i*-zr+ E (**)0"*i) (3-33) 

i+j=m UJ ^i u -Lj —i-\-j=m j —i—j=m 

j (3) _ V 7 (1) 7 (1) ■ 

<J m — • °i °j °k 

i+j+k=m 

<9 3 d 2 
= V I 3 V ix 



l+J — k=m dxidxjdxk _ i+J+k=m dxjdx k 

d 

+3 E ( ix i)ti x j)n E {ixi)(jxj)(kx k ) (3.34) 

. r-f, oxh ■ , 



where i,j, k = 1,2, 3, .... 
By setting 



1 

n 2 
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we find 

J {2 \ = 2t *i + t\ + 4 ^ ( n+ I)t n+2 _ (3.36) 

n=0i,l, " 



^- 2 2=^ + 2 E (2n + l)t n+1 J- (3.37) 

n=0,i,l- ^ 

^=4 E + ( 3 - 38 ) 

n=0,i,l,.. 

4 2) = 4 E + ( 3 - 39 ) 

nil ^ 



From (I3.34p . we have 



a: 

jf) = StlU+3 V (2j + l)(2fc + l)t 



+ s E + ( 3 - 4 °) 

Then, the first equation for the Virasoro constraints is expressed by 
The second equation becomes 

2 

The third equation is expressed by 

3 wrr,-\ e + kt °H ~ \ J °' 2) ' i J -> m )» = (3 - 43) 

The differential operator appears only for the equation of order A^ n 
(n=l,2,3,...). This is similar to the p-spin generalized Kontsevich model 
without logarithmic term, where spin equations are described by J^f 1 and 
the spin non-zero equation of motion is described by jffl [33J. 

If we denote the differential operator \ J2m^ as L m : 

L n = \j% (3.44) 
those L n have the commutation relations 

[L n , L m \ = {n- m)L n+m (3.45) 
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4 Integral equation for the Airy matrix model 



For the unitary matrix model the large N limit may be solved by a Riemann- 
Hilbert integral equation [IT] . We apply here the same technique to the 
Kontsevich model with a logarithmic potential. 

When k=0 (the Kontsevich model), the equation of motion reduces to a 
simpler second order equation. Let us first consider the k = case as an 
exercise, following [llj. 



dX 2 . ^ X c - X d \dX c dX _ 
Changing to the free energy W 

Z = e PW (4.2) 
(the original Gaussian matrices were NxN, but the dual matrices are PxP), 

Introducing the density of eigenvalues 

p(*) = 4 (4.4) 

we consider W as a functional of p from which one obtains W a as 

Wa = w(x)\ x=Xa (4.5) 

with 

w f x ) - 1 d 5W (4 6) 

P dx Sp(x) 

The second derivative in (14. ip leads to two terms, but in the large P-limit 
the leading one is simply w(x) 2 , leading to the integral equation 

9/ x f b , /x w(x) — w(y) 

w 2 (x)+ dypiy)^^ — = x 4.7 

J a x — y 

We define / and F as 

f(z)=[ b dx-^- (4.8) 

J a Z — X 

F(z)= f dx p{x)w{x) . (4.9) 

J a Z — X 

Inside the cut z G [a, b], 

ReF(z) = w 2 (z) + w(x)Ref(z) - z (4.10) 
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and for z G [—00, 00], 

ImF(z) = w(z)Imf(z) (4.11) 

We make the ansatz 

F(z) =w 2 (z)+w(z)f(z) -z (4.12) 

leading to 

Imw(Re/ + 2Rew) = z G [-00, 00], (4.13) 

(Iimu) (Imw + Imf) =0 z G [a, b] (4.14) 
From (I4.12p . in the z — > 00 limit, we find F ~ 1/z, f(z) ~ 1/z, and 

w ( z ) = v ^--L + 0(^ 3/2 ) (4.15) 
Since Ymw 7^ for 2 G [—00, — c], we have from ( 14.131) . 

Rew = ~2 Re f ( 4 - 16 ) 

This is equivalent to 

lm(w(z)y/z + c) = -^f(z)y/-z - c (z G [-00, -c]) (4.17) 
Then by dispersion relation, we get 

1 r, f(y)V=y^ 



1 r 

w(z)y/z + c = --J dy- 



z-y 



± r d y[ b d X P {x) ^~\ (4.18) 
Z7T J-oo J a (y — x)(z — y) 



Noting that 



dy , ^T~ C =Tdt- ^ 



00 (y — x)(z — y) Jc (t + x)(t + z) 

71 



y r z + c+ y/x + c 



(4.19) 



Adding the integral constant z+~, which is determined from the asymptotic 
behavior of (I4.15p . we get 

c 1 f b p{y) 



w{z)Vz + c = z + --- dy^=^—= (4.20) 
2 2 J a a/2 + c + y/y + c 

and the parameter c is determined from the condition that there is no pole 
at z — — c, 

rb dy^L (4.21) 
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Thus we obtain 

w(z) = ^fz~T~c + - [ b dy , , P( ^^ ; (4.22) 

2 J a {sjz + c + y/y + c)yjy + c 

This function w(x) is indeed a solution of the integral equation (14. 7p . The 
square of the second term, the part of the integration in (I4.22p , cancels with 
the second term of ( 14. 7ft . and the parameter c is given by (14.2ip . 

By further integration over p(x), we find the free energy W, 

2 fb 3 rb 

W = - / dzp(z)(z + c)2 — c / dzp(z)\/z + c 

O J a J a 

- t dz f dyp(z)p(y)\og(Vz~T^ + v 7 ?^) - ^c 3 (4.23) 
where we used 

dW d 5W{p) 

3K = S W = ™ (2) <424> 
The parameter c satisfies the saddle point equation for tu in (I4.23p . 

dc 4 Y v Ad + c 

and this is consistent with (14.2 ip . 

The large P limit of the free energy F in ( I3.30P is obtained by the scaling 
Aj ~ P 2 / 3 and by taking each sum as order P. The order of t n becomes 

*» = E-^~°( pl(1 ~ ft) ) ( 4 - 26 ) 

The expansion of c in (14 .251) is obtained by the recursive solution with 
the definition of t n (t n = 2+n ), 

c = -E 1 



\Ai + c 

= -t - ^oti - ^2 - 7^0 + 0(A" 5 ) (4.27) 
z o 4 

From this equation, c has to be negative, and A > — c. Therefore, we have 
only one expansion, the large A expansion. 

The free energy W = logZ in ( 14 . 2 3 [) is divided into four terms. We expand 
each term for small c (c is a constant expressed by t n ), 

W 1 = §£(A, + c)i 

15 



2 3 i c 2 1 c 3 1 

3 « i t A, 2 24 t A, 2 



W 2 = -cE(A, + c)^ 



i r 2 i 3 i 

-cEAf-^E J r + ^E J 3+- 



2 A 2 - , A/ 



w 3 = -^E lo g(A + c + \/^+ c 

= -^EMv / ^+v / ^)-iE 7r |TT + |E4T + --- 

W 4 = -- (4.28) 

Inserting the expression of c, we obtain W, which is the sum of these four 
terms, 

W = I E A? - \ E + \/^) + + ol*o*i + O(i) (4.29) 

i ij A 2 

The first two terms are log-Zo i n (13. 5p and remainings are consistent with the 
genus zero part of F in (I3.30p . Up to order 1/A 3 , only t 3 , and tjjti are genus 
zero terms. 

The free energy F is of order P 2 in the large P limit. From (I3.30p . we 
find in the large P limit, 

d 2 F 11 3 1 

U = W = 2 to + 4 totl + lQ t2 ° h + 8*°*' + " ' (4 ' 30) 

Thus we find 

d 2 F 

c = -2u = -2- w (4.31) 

Therefore, we understand that c is the specific heat for the free energy 
F, when we interprete to temperature. 

We now consider the Kontsevich model with a logarithmic term (k ^ 0). 
The equations of motion in (13 ,4p are expressed as equations for W and W a , 
where 

In the large P limit, the first and the second terms in (I4.32p are negligible. 
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From f )3.4p . we express it as 

■\, \ , \ f w(x) — w(u) , . . . , , 

w [x) — xw(x) + 2 — -p[u)p{v)dudv 

J (x — u)(u — V) 

-(1 + -)+ f du^-(2w 2 (x) - w(x)w(u) - w 2 (u)) = (4.33) 
P J x — u\ J 

From f )4.33p . we find in the large x limit, 

w(x) ~ y/Z - ±(1 + ±) (4.34) 

This is a generalization of ( 14.1 5ft for k ^ 0. 

The equation of f 14 . 3 3 1) is a cubic equation. If w(x) has a solution similar 
to (14. 22 p . the tri-linear terms of p has to be cancelled in this cubic equation 
of w(x). First we check that whether the solution of (14.221) satisfies (14.331) 
when k = 0. We denote 



l x = y/x + c, l y = yjy + c, l z = V ' z + c, / s = \fs + c (4.35) 
The solution for /c = is 

^^4/%^ < 436 > 

We express the cubic equation of (14.331) in terms of these l x ,l y ,l z , l s by 

Ji + I 2 + / 3 + h + / 5 = (4.37) 

3 12 f J„. 



A = w \ x ) = ll + -llJ dy 
4 /" 



(^x "I - ^y^y 



(l x -\- ly){l x ~\~ lz)lylz 

+ - [ dydzds P(y)p( z )p( s ) ( 4 

I 2 = -xw(x) (4.39) 
2 / dydz K \ yy> p{y)p{z) 

J {x-y){y-z) 

2 / dydz- x — * -p{y)p{z) 

+ /.frfa* . ^M') ( 1 - 1 ) (4.40) 
J (x - y)(y - z)l s l x + l a ly + ls 
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h = -l 



(4.41) 



dy-^- ( 2w 2 (x) - w(x)w(y) - w 2 (y)) (4.42) 
x — y V / 



Up to the first order of p, by adding the contribution / of Ji, I2, h, h, h, we 
have 



AI 



p{y) 



U dy ] 

-\ J d v dz 7p^n ( 4 - 43 ) 



p{y)p{z) 

(l x ~\~ lyjlylz 



where we have used the expression of c given by (14.2 lj) . The summation of 
above AI and the second order of p in I\ and J5 becomes 

AI + I 1 + h = 2f dydz .. P }f$ z) (4.44) 
J {lx + ly){z-y) 

and this is cancelled by the contribution of I 3 . Thus the contribution up to 
the second order is cancelled. The terms of the third order of p come from 
I\, I3 and I5. There are triple integrals over y, z and s. We symmetrize the 
integrals over these three variables. Before making the symmetrizations, we 
note that ()()() 

8 l ^ (l x + l y )(l x + l z )(l x + l s )lylzh ^ ^ 

(ii has symmetric form ). 

h = - I dydzds P{y)p{z)p{s)(l* + ly + h + ls) (4 m) 

2 J (l X + l S ) (l X + ly) {ly + l S )(l X + Q{l Z +l s )(l y +l z )l s ^ ' ^ 

h = -- f dydzds P(y)p(z)p(s)(l x + l a + 2l y + 2l z ) 

4 7 (l x + l y )(l x + l z )(l y + l z ) (l x + l s )(l y + i s )i z i s v • ' 

After symmetrization, these three terms cancel completely (/1 + /3 + J5 = 0). 

Thus, we see that the equation (I4.33j) is satisfied by the solution of w(x) 
in tKTZty when k = 0. 

Since the parameter k appears only in terms of order zero of p in (I4.33p . 
it is easily understood that there is a straightforward solution for k 7^ 
case, based on the above analysis. Since we have seen the solution of (14.221) 
satisfies the integral equation of (14.331) . we consider the solution (14.221) for 
k 7^ more carefully, specially the condition for c. Namely, we use the same 
solution w(x) as before 



w(x) 



Vx~+~c + l -( b dy___4?L_ == (4.48 

2Ja y (y/xTc+ ^Y+C)V¥+C 
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where we consider that the parameter c is now a function of x, c = c(x). We 
replace all parameters c by c(x) in (14.48j) . Since x is fixed in the integral 
equation, this change from a constant to x-dependence of c does not make 
any difference. 

Up to first order in p, by putting this w(x) into the part of first order in 
p of f)4.33p . we have 



w 3 (x) — xw(x) — (1 + — ) + / dy ( 2w 2 (x) — w(x)w(y) — w 2 (y) 

P J x — y \ 



= (x + c) 2 — Xy/x + c — (1 + 

3 . f b p(y) 
+ (x + -c" 1 1 



2 Ja (y/x + C + y/y + c)y/y + C 

'• p{y) 



+ [ dy-^- ( 2(x + c)-Vx + c^yTc - (y + c)) 
Ja x — y\ J 

k i — i— r b , p(y) , c r b p(y) 



c^x~T~c -- + yfx~Vc j dy-^= + £ [ 

P J a \/y + C 2 Ja 



(4.49) 



where we have used that the integral of p is one, 



p(y)dy = 1 (4.50) 

J a 

We now put 

c= - t dy-?ML + h (4.51) 
J a ^yTZ 

and r.h.s of (I4.49P becomes 

k h f b p(y) 



r 



h.s = h\/x + c - — + - j 

P 2 J a 



P 2J a (Vx + c+ y/yTc)y/y+~c 



1 f dz P ^ fb P ^ 



2J y/z + cJa (y/x + c + y/y + c) y/y + c 

(4.52) 

The last term is transfered to the part of second order in p, and we have 
seen that the second and the third order terms of p are cancelled completely. 
Therefore, we find that 

1 fb p(y) 



(y/x + c+y/y + c)y/yTcJ P 



k 

hw(x) - — = (4.53) 
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Thus we have the coupled equations 
c = - f dy 

J a 



p(y) k 1 



yjy + c Pw(x) 
w (x) = Vx~T~c +\f dy pi -fL_ (4.54) 

2 Ja (y/X + C + y/y + Cjy/y + C 

This w(x) satisfies the integral equation (I4.33P and in the large x limit, it 
satisfies the asymptotic behavior 

1 k 

w(x) ~ yfiH + — (1 + — ) (i4oo) (4.55) 
2x j 

We obtain the expansion of w(x) for large x from the coupled equations 
in f)4.54p . By integration over p(x) we obtain the free energy W for large A 
(■'• A). 

+ ^o + ^i + ^ 2 t 1 + --- (4.56) 

which is consistent with logZo in ( 13.51) and the genus zero part of F in (13.301) . 

When k is sufficient large, we have a solution in which c is positive in 
(I4.54p . In this case, we obtain an expansion for small A. 



5 Intersection numbers in a replica limit 

In the case of one matrix model, we have used a duality relation between 
the Kontsevich model and the Gaussian random matrix model at a critical 
edge point [U EJ |3]. More precisely the Fourier transform of the n-point 
correlation function U(si, s n ) becomes the generating function of the in- 
tersection numbers with n-marked points. This n-point correlation function 
U(s\, ...,s n ) has a Cauchy integral representation, which is equivalent to the 
integral of the first Chern class over the moduli space -M 9) „. 

We have shown in section 2, that there exists a similar duality relation 
between the partition function of the Kontsevich-Penner model (II. ip and 
the correlations for the Gaussian distribution ( I2.5p . We want to discuss the 
origin of the terms (half-integer) in this section. 

In the expansion of the free energy F in (13 .30 p . the number of times of 
appearance of tj is the number of marked points according to the definition of 
the intersection number in (12. ip : In ( 12. ip . n is the number of marked points. 
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We investigate first the case of one marked point (i.e. a single tj). This 
case is obtained from the replica limit for the matrix B in ( II. ip . namely the 
limit in which its size P goes to zero. We first make a shift B — > B + A 2 to 
eliminate the linear term trBA, 

Z=[ ^ e -|trB 3 -trB 2 A5+fetrlog(A3+B) /q-q 
J PxP 

where B is a P x P Hermitian matrix. The replica limit for B means that 
we take P — > limit, selecting thereby the contribution for one marked 
point [21 H]. In this replica limit, the number of eigenvalues Aj also goes to 
zero, and it becomes not necessary to distinguish them. Therefore, we set 
A = A • I. Then, the calculation becomes considerably easier. First, we make 
the rescaling B to B/X 1 ^. Then the partition function Z in (15.11) becomes 

Z = J djBe -^-itr B 3_ trB 2 +fctrlog(A | +B) ^ ^ 

We first neglect the cubic vertex tr£? 3 . Then, one recovers for Z the 
same model when p = 1 in (jl.2p . as was discussed in references [131 E3]- We 
consider here this same model by the duality plus replica method. 

Z = J dBe' trB2+ktTlog ^ +B) 

= 2-£ / dB[det(Al + 2- 1 iB)} k e-^ trB2 (5.3) 

JPxP 

From the duality theorem for characteristic polynomials |2j, the above 
expression has a dual form, which is an integral over a k x k Hermitian 
matrix M : 

Z=( dM[det(M + At)] p e -5 trM2 (5.4) 

Jkxk 

In the limit P — > 0, we have (A = v^A^) 

lim = / dMtr^^e-^ trM2 (5.5) 
p^o ^xfe A + M 

which is the one-particle Green function G(X) for the Gaussian random ma- 
trix, and the expansion of the inverse of A is easily obtained as a moment of 
M in terms of polynomials of k. 

By integration of G(X) about A, we obtain Z, 

Z = klogX + y < tiM 2j > (5.6) 

U(2j)X^ 
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The Gaussian average < trM 2 - 7 > is easily evaluated from the integral 
representation of U(s) [2], 

U(s) = < tre sM > 

= ! e W*!-(i + £)* e « (5.7) 

s J 2m y u' y ' 

This integral runs over the contour centered at u — 0, and it yields 

,,, , . k 2 2 2k 3 + k 4 fc 4 + 2fc 2 6 2fc 5 + 10P + 3k 8 . . 

^U^+T' + ^r s + -[44~ s + — 5760 — s + "-- (5 - 8) 

From these expressions, we obtain 

< trM 2 >= k 2 , < trM 4 >= 2k 3 + k, < trM 6 >= 5A; 4 + 10A; 2 , 

< trM 8 >= Uk 5 + 70k 3 + 21k, ■ ■ ■ (5.9) 

From (15. 6p , we express Z of (15.61) in terms of tj (tj = — — r), by noting 
that s = j = — -fhr, 



A; 2 1 1 
Z = k\og\ + —t x + —(k + 2k 3 )ts + —(5k 4 + 10k 2 )U + ■■■ (5.10) 
4 16 2 48 

The term \k 2 t\ coincides with the term in F of fl3.30[) . 

The tri-valent term, which we have neglected, couples to the logarithmic 
term, and also make a contribution as polynomials in k. The exponent 
exp(— |A _ 2tri? 3 ) is expanded and it gives the contribution in the replica 
limit P — > 0. By the formula of the replica limit, we have nonvanishing 
average of < n« trl? di > . This formula is [2] 

Iiml7(si,...,sj) = lim 4 < tre SlB ---tre SiB > 
p^o v ' ' 11 p^o P 

= llpsinh^ (5.11) 

° 3=1 2 

where a — s% + ■ ■ - + S/. This provides a generating function for < FL tri? di >. 
From this formula, for instance, we have 

lim — < tr^tr^ 3 >= 3, lim — < trB^r^trE 2 >= 18, ■ • • (5.12) 

P->0 P P-s>0 P y ' 

Using these values of averages, we are in position to compute the co- 
efficients of the t n terms. We consider the term ts in (13.301) . which has a 
coefficient \(k + k 3 ). The partition function Z is 

Z = j dBe~^ tlB3 ~ tlB2A ^ +ktTlos ( 1+A ^ B ^ (5.13) 
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We rescale B — > 2 s A J 5 and A 2 — > 2 3 A 2 . We have 



/ rfSe -^ trB3 ^ trB2+fetr(1+A " lB) (5.14) 



Expanding then the logarithmic term and exp(— 6A 3y 2 tr-B 3 ) term, we find the 
contributions to order A from 6 terms. These terms are evaluated by the 
replica formula (P — > limit) of (15.111) . 

/ .\ 1 k „/i k 

1 1 o/c,l*o„o A/ 



P 6A* 3 A* 6A4 

WF<^'W-|(> rB2 >=8^ 

Adding these (i)~(vi) terms, and noting that we have made a scaling of A, 
we obtain as expected the result |(fc 3 + A;) A -3 = |(/c 3 + k)ts for f)3.30p . 

We now evaluate the coefficients of tj from the Fourier transform of 
one-point correlation function U(s). We have shown that the intersection 
numbers of n-marked point are obatined from the Fourier transforms of n- 
point correlation function U(si, s n ) [Sill]- The Kontsevich-Penner model 
involves two interaction terms tr5 3 and trlog.8. For the application of 
U(si, s n ) to this Kontsevich-Penner model, we have to extend the pre- 
vious duality expression. 

Let us return to the duality relation for the one-matrix model, before 
extending it ; in this one-matrix case the duality reads 

k N 

< [] det(A Q - M) >m,a=< II det ( a i - iB ) >b,a (5.16) 

a=l j=l 

in which the l.h.s. consists of Gaussian average in an external matrix source 
A for Hermitian N x N matrices ; the r.h.s. is also a Gaussian average for 
Hermitian k x k matrices in an external source A (whose eigenvalues are the 
X a of the l.h.s. ; the aj are the eigenvalues of A). We also know in closed 
form the Fourier transform of the n-point correlation function, 

n 

U( Sl ,...,s n ) =< Y[tre s ' M > M>A 



23 



= fjr dUi e NX«t>i+kNE* ft TT(1 + — )det (5.17) 

J ZZ71 f-J^ J = 1 U\ — Clj Ui — Uj + Si 

Let us consider the one-point function 

= - / -^exp(- £ + s) m - u m ] + iV MS + -Ns 2 ) (5.18) 



1 * 



where 

c m = — V-^— (5.19) 
The r.h.s. of the duality formula (15. 16ft is also expanded in powers of B, 
< JJdet^-zfl) >= j dBexp(-Y,c m (iB) m -^NtrB 2 + NtTBA) (5.20) 

From this representation we have investigated the (p,l)-model (the (2,1) 
corresponds to Kontsevich model), which is obtained by specifying appro- 
priately the a,j [3]. We consider here a more general situation, in view of 
encompassing the Kontsevich-Penner model. The (p,q)-model is defined by 

Z = J dBexp(~c p+1 tiB p+1 - c q+1 tiB q+1 + txBK) (5.21) 

This is obtained by imposing the following conditions to the aj, 

1 * 1 N 

2 h ~ 2 



1 n 



E TTw = ( 5 - 22 ) 



where m = 3, 4, q and + These conditions should be understood 

as holding in the large N limit. 

These conditions may also be applied to U(s) in ( I5.18p . Then we have 

For the application to the Kontsevich-Penner model, we have to take the 
limit p — > — 1 and 0—7-2. Since c p+ i = -^7 J] -, — k+r. we have 

f P+L \ a j) v 

Z = f dBexp(-c 3 tiB 3 - iVtrlogE + trSA) (5.24) 

Jkxk 
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U( S ) = 11 *i e - C3 (( U+S )3- U 3 ) _ mrlog( ^ ) (5 _ 25) 



and 

U(s) = , 

s J 2m 

For the n-point correlation functions, we have similarly 

U( Sl ,...,S n ) = /n^i e -^ C 3(K^) 3 -f)-iVE i trlog(^) det 1 

J fJl 2«7T Ui - Uj + Si 

(5.26) 

The intersection numbers for one- marked point, which are obtained from 
Z as the coefficients of the linear terms in the tj 

Z = < Tj > tj + higher degree, (5-27) 

are derived from U(s) as coefficients of the expansion in powers of s . 

Apart from notations, in which we have to interchange N to k, and further 
k —k, and chose C3 = |, Z is then identical to the Kontsevich-Penner 
model in (11. ip . We then have 

17(*) = e — I ^e-^C^-Y (5.28) 

S J 2m u 

We now expand U(s) in powers of s andl k. To that purpose we shift u = 
v — |s, and v = sz/2, and obtain 

This gives an expansion in powers of k, from 

i^r? = l + Mog(^±|) + ^ 2 [log(^±|)] 2 
z — 1 z — 1 2 z — 1 

+ ^![log(^±i)] 3 + 0(fc 4 ) (5.30) 
o z — 1 

The first order of leads to the integral 

/°° 1 1 o 7 A- 1 7T q 1 ^ 

cfee~* s * log(^-) = 2i(-)3er/(^si) (5.31) 
-00 z — 1 s 2 

where erf(x) is the error function. For small s and small k, we obtain 

U(s) = k - Iks 3 + 0(ks 6 ) (5.32) 
6 

The integrals for odd powers of the logarithm may be computed analytically. 
For instance, from ( 15. 29ft . in the the small s and k expansion, we obtain a 
term of order of s 3 k 3 , 



3 k 3 t°° , .1 z 2 .l n ,z + l N1 

L dz( T2 + w ios( — )] 



3 



2 



lms 3 k 3 (5.33) 
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Therefore we recover from 05.321) and (15.331) . the coefficient of s 3 as \{k + 
k 3 )s 3 , which yeilds the expected + k 3 )ts term in (13. 30 h . 

We note that the odd powers of k may be obtained systematically by 
computing a residue at z — 1 in the contour integral (I5.29p . For the even 
powers of k, we have to integrate the logarithmic terms. 

The replica formula (15.111) has been derived with a single logarithmic 
integral, which has a cut for — s < u < [2J. This replica formula corresponds 
to the k — > limit of the(p,q)-model with q = 1 and p = — 1, which leads to 
the intersection numbers with one- marked point. In an appendix, we present 
the two-point function U(si,s 2 ) for q = 1 and p = — 1 as a polynomial in 
k, consistent with the replica formula for k — > 0. From this example, we 
understand that the contours for the n-point function should encompass all 
poles. We now consider the two-point case for the Kontsevich-Penner model 
with two marked points : 



N i (s 3 +s3) r d Ul du 2 e m s 1+Sl u 1+ u 2 s 2+U2 s 2 u +s u +s 
U(s u s 2 ) = e ^ +s *> f — — — — 

J (2'ITT) 2 [Ui - U 2 + Si)[Ul - U 2 + S 2 ) Ui u 2 

(5.34) 

where the contours for u 2 are such that one sums over the three contributions 
from the poles at u 2 = 0, u 2 = u\ + s\, u 2 = u\ + s 2 and then the contour for 
u\ circles around the origin. This yields 

k 2 

U(s 1 ,s 2 ) = k 2 (s\s 2 + sisl) + — (sls 2 + s lS 5 2 ) 

+ 4 + T2 k2){sisl + sl4) + ( T + l k2)sl4 + °(^ 5 - 35 ) 

Changing from s m to t m _i, we obtain the terms with two marked points of 
(ESQ]). 

Thus we find that the terms with n-marked point for the Kontsevich- 
Penner model (11.11) are expressed explicitly by the integral formula of U(s\, s n ) 
of (I5.26p . The Kontsevich-Penner model is the p = —l,q = 2 of the (p,q)- 
model of f 15 .211) . Following the same techniques as above we find similarly ex- 
plicit integral representation of U(s\, s n ) for the (p,q)-model with p — — 1 
and arbitary q. 



6 Discussions 



In this paper, we have extended the analysis of our previous paper j3] to 
the Kontsevich-Penner model. We have derived the Virasoro constraints for 
this model, and we have obtained the large N solution of the corresponding 
integral equation. The occurence of parameters t n , with half-integer n, is due 
to the logarithmic potential. Using the correlation functions of the Gaussian 
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two-matrix model, in a source, with one set of eigenvalues near an edge, and 
the other one in the bulk of the spectrum, provides the Kontsevich-Penner 
model. 

We have used an explicit integral representation for U(si, s n ) (15. 26 p . 
which gives then fc-dependent coefficients of the free energy F of the Kontsevich- 
Penner model. This model turns out to be the special limit p — — 1 and q = 2 
of a (p, g)-model . The integral representation for U(si, s n ) is valid also 
for general q with p = —1. The details for such cases is left to future work. 

In string theory, the c = 1 matrix model has attracted considerable in- 
terest, renewed recently from the D-brane point of view. The tachyon plays 
a central role in the c = 1 matrix model. In the present study, the partition 
function for the Kontsevich-Penner model ( II. ip . is derived from a time de- 
pendent Gaussian matrix model, and the role of the and fc-dependence are 
clearly understood from the correlation functions of the two-matrix model. 
Thus, it may shed a light on the c = 1 string theory, FZZT-brane, etc [401 . 
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Appendix A: formula for p-th derivatives 



The matrix A has eigenvalues Ai, A2, .... and corresponding ortho normal eigen- 
functions \(fi a )- We now consider a perturbation dA, 



(A + dA) (J0 > + \d(fi > ) = (A + dX) >+\d<f)>^ 

and from this equation, we obtain at first order 

(A - X a )\d(fi a > +(dA - dX a )\(fi a >= 
Multiplying < (fi a \ from the left side, it becomes 

< (p a \dA\4> a >= dX a 
In an arbitrary fixed orthonormal basis \b >, it becomes 

dX a =< (fi a \b >< b\dA\c >< c\(p a > 
Therefore, we obtain the first important formula, 

dX n 



OA 



- =< <f) a \b >< c\(j) a > 



be 



< <f) b \d<f) a > 



X„ — Xh 



< 6 h \dA\6 n > 



(a.i; 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



Note that < (p a \b >= U a b, where U is a unitary matrix. From ( 1A.2j) . multi- 
plying by < (fib | (b 7^ a) the left hand side, 



(A.6) 



Therefore, we have 



\ d( Pa >= E 



X a — Xb 



b b >< 4> b \dA\0 a > 



from which follows the second formula, 



d < b\cfi a > 



dA 



cd 



E 



The conjugate of this formula is 



d < (fi a \b > 



OA 



dc 



E 



1 



A a — A/ 



(A.7) 



< b\cfif >< (fi f \c >< d\(fi a > (A.I 



< (fi f \b >< c\(fi f >< (fi a \d> (A.9) 



28 



By the chain rule, we obtain the first derivative, 

dZ d\ c dZ 

dA ab dA ab d\ c 

= < b\<f) c >< <p c \a > i^—j (A.10) 

The formula for the second derivative is obtained by the use of ( 1A.5|) and 
([XI). 

Z = (°A (°A Z 



KdAVab \dAJad\dAj 



db 



<} f <b\<t> c x<t> c \d> f|^) ) (All) 



dA ad \ " C| ' V^A 

Noting that 



< d> r \d > 



)flAZ <6| ^ > 



<9A 



and 



; < <Pc\d > (— ) ^ - _ A < 6|0/ >< 0/|a >< d|0 c > 

eye / 

:?< ^ c>< ^i a> (^)_i_ (A . 12) 



< **- > (wi)<L <Md> 



=< b\(p c >< (p c \a > . . ^ 

\OA c / d A c — Ad 

we obtain 



The third order differentiation is obtained by repeating the same proce- 
dure. We write T c by 

Tc = ^7 + £ ~ W d ) (A - 15) 
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(JL) = (±) (£.) 

iMl)ik( <b ^ >T < <Ma> 



'■pa 

or 

= < b\(f) c >< 4> c \p >< a\(j) c >< (j) c \a > 



c 



+ < (f> c \p >< a\(p c > T c < <j) c \a > 



d\ c 
d < b\<j) c > 



+ < c |p >< a|0 c >< 6|0 C > r /< ^l a> (A .16) 



Therefore, we obtain 
dWJab 



< b\<Pc > + E y^r^ - r ^)) < > ( A - 17 ) 



Using the identity, 



1 1 1 / A \ 

+ 7^ \~~\7\ + \~~\7\ TT = ( A ^) 



(A c — A d )(A c — A e ) (Ad — A c )(A d — A e ) (A e — A c )(A e — \d) 

we obtain the expression in terms of eigenvalues 

d 3 \ d 3 
dX 3 )^ ~ 9A? 

+ ^ C A C -A/9A: _ W ( 9A: + 0V " ^ (A c - A d ) 2 { d\ c " w d } 

j^ c ^ (A c — A e )(A e — A d ) ^9A C 9A e ^ ^ ^ 
If we write 

r?> = ± (A.20) 

we have 

( M^) ab =< 6|0C > r ^ < 0c|a > (A ' 21) 

r? } = + £ r ^ T (r« - if) (A.22) 

Repeating this procedure, we obtain 

r&H-D = 9 r w + E - J_(rw - if) (A.23) 
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and 

(aA^i) a6 =< ^ > Ti c +1) < &|a > (A.24) 
Appendix B: Relation to unitary matrix model 



We will show that the unitary matrix model with external source [TT] is 
equivalent to the higher Airy matrix model with a logarithmic potential for 

p = -2 0Eg. 

The unitary matrix model (Brezin-Gross model [TT]) is 



Z = J dUe^ UAt+u ' A ^ (B.l) 



where U is a unitary matrix and A is an arbitrary complex matrix. From 
the unitarity condition, UW = 1, we have 



d 2 

-Z = I-Z (B.2) 



Introducing A as 
we find 



DAW A 

A = AA* {A l3 A) k = A lk ) (B.3) 



d 2 d dA qs d 



Thus 



3A\ 3 dA ]k ~ dA% dA jk dA qs 
d d 
dA ks ms dA mi 
d 2 d 



x & +n m) z=i - z (b - 5) 



The equations of motion for the p = —2 case follows from 



dM JL e trM-i+trMA+ktrlo S M = g ( B . 6 ) 

Taking two derivatives with respect to A, this leads to 

l + 2N—+A—-+k—)Z = (B.7) 
OA OA 2 dAJ v ' 

If we now take k = —N, ( \B.7\i this equation becomes identical to (IB . 5 P . 
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Thus we find that unitary matrix model with an external source is similar 
to the Kontsevich model with a logarithmic potential. There is a phase 
transition in this unitary matrix model with a critical point at 



1 A 1 



tr^= = V^= = 2. (B.i 
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